Negative Index Lens Aberrations 
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We examine the Seidel aberrations of thin spherical lenses composed of media with refractive 
index not restricted to be positive. We find that consideration of this expanded parameter space 
allows reduction or elimination of more aberrations than is possible with only positive index media. 
In particular we find that spherical lenses possessing real aplanatic focal points are possible only 
with negative index. We perform ray tracing, using custom code that relies only on Maxwell's 
equations and conservation of energy, that confirms the results of the aberration calculations. 



In 1968, V. G. Veselago proposed the idea that a ma- 
terial could have a negative index of refraction, and de- 
scribed how this would impact many basic electromag- 
netic phenomena 1]. In recent years, there has been great 
interest in this subject due to experimental demonstra- 
tions of negative index artificial materials^, and the in- 
troduction of the perfect lens concept fs!]. A perfect lens 
is a flat slab of index minus one, which can focus images 
with resolution exceeding that possible with positive in- 
dex optics. 

Recently, focusing by a negative index medium 
with curved surfaces has been experimentally 
demonstrated^. Traditional spherical profile lenses 
composed of negative index media have several advan- 
tages over their positive index counterparts: they are 
more compact, they can be perfectly matched to free 
space, and here we demonstrate that they can also have 
superior focusing performance. 

The monochromatic imaging quality of a lens can be 
characterized by the five Seidel aberrations: spherical, 
coma, astigmatism, field curvature and distortion. These 
well known corrections to the simple Gaussian optical 
formulas are calculated from a fourth order expansion of 
the deviation of a wave front from spherical. (A spher- 
ical wave front converges to an ideal point focus in ray 
optics). The coefficients in this expansion quantify the 
non-ideal focusing properties of an optical element for a 
given object and image position^. We find that there 
is an asymmetry of several of the Seidel aberrations with 
respect to index about zero. Considering that an inter- 
face with a relative index of +1 is inert and one of rela- 
tive index -1 is strongly refractive, this asymmetry is not 
surprising. However, our conclusion that the asymmetry 
can yield superior focusing properties for negative index 
lenses is not obvious. 

We note that negative index media are necessarily 
frequency dispersive, which implies increased chromatic 
aberration and reduced bandwidth. However, diffractive 
optics, which possess a similar limitation, have found util- 
ity in narrow band applications]^. 

To confirm the analytical aberration results, we devel- 
oped custom ray tracing code that does not rely on the 
sign of the index to determine the path of the ray, but 
relies only on the permittivity, £, the permeability, ^, 
Maxwell's equations and conservation of energy. 

Between interfaces, in homogenous media, the ray 



propagates in a straight line following the direction of 
the Poynting vector. Refraction across an interface, from 
a region labeled 1 into a region labeled 2, is handled as 
follows. Wave solutions are sought that satisfy the dis- 
persion relation (obtained from Maxwell's equations) in 
region 2, 



— k2 • k2 = e2Ai2, 



(1) 



where k2 is the wave vector in region 2. The solutions 
must also satisfy a boundary match to the incident wave, 
requiring 



n X (k2 - ki) 



0, 



(2) 



where n is the unit normal to the interface. The outgoing, 
refracted, wave must carry energy away from the surface 
if the incident wave carried energy in, 



(P2 • n) (Pi • n) > 0, 



(3) 



where P = ^ Re (E x H* ) is the time averaged Poynt- 
ing vector. Finally, the wave must not be exponentially 
growing or decaying, Im(k2) = 0, since the media are 
assumed passive and lossless. If a solution exists that 
satisfies all the above criteria, the ray is continued with 
the new found wave vector and Poynting vector. Further- 
more, since we consider only isotropic media the solution 
will be unique. 

We find that the form of the expressions for the Seidel 
aberrations of thin spherical lenses found in the optics lit- 
erature are unchanged by the consideration of negative 
index media. We reached this conclusion by re-deriving 
these expressions, from first principles, using only the def- 
inition of optical path length and Fermat's Principle. We 
interpret the optical path length, OPL = J^n(s)ds, to 
be the phase change (in units of free space wavelength) 
that a wave would undergo along the path C, if C is 
oriented parallel to the Poynting vector. The optical 
path may have contributions that are negative where the 
Poynting vector and the wave vector are antiparallel, i.e. 
where the index is negative. These aberration formulae 
are further corroborated by agreement with the results 
of our ray tracing. 

The wave aberration, AOPL, is the difference in op- 
tical path length of a general ray and a reference ray. 



2 



general ray 




FIG. 1: Construction used for aberration calculation. The 
aperture stop, labeled AS is at the plane of the thin lens 
(though lens shown is thick). The Gaussian image plane is 
labeled IP. The aperture stop coordinate vector, r, and the 
image plane coordinate vector, h, are not necessarily parallel 
as shown. 



where the reference ray passes through the optic axis in 
the aperture stop and the general ray is parameterized by 
its coordinate in the aperture stop, r, and its coordinate 
in the image plane, h (Fig. To be in the Gaussian 
optic limit, where spherical interfaces yield perfect imag- 
ing, r and h must be near zero. A series expansion of the 
wave aberration in these parameters 



AOPL= G„„(r.r)'(r.h)"(h.h) 



(4) 



yields corrections to Gaussian optics of any desired order. 
The lowest order corrections for a thin spherical lens with 
aperture stop in the plane of the lens are given by 
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(5b) 
(5c) 

(5d) 
(5e) 



These coefhcients are the Seidel aberrations: spherical, 
coma, astigmatism, field curvature and distortion respec- 
tively. Also appearing in these expressions are p, the po- 
sition factor, and q, the shape factor, where we follow the 
definitions of Mahajan 5]. The position factor is given by 



2/' 
S' 



(6) 



where /' is the focal length referred to the image side 
and S' is the image position. Through the thin spherical 



lens imaging equation. 
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(7) 



where S is the object position and i?i and i?2 are the lens 
radii of curvature, the position factor is directly related 
to the magnification. 
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The shape factor is given by 



i?2 + Ri 
i?2 — ^1 



(8) 



(9) 



A lens with a shape factor of is symmetric, and ±1 is a 
piano-curved lens. Using the shape and position factor, 
all thin spherical lens configurations are described. 

We will first examine the very important case of a 
source object at infinite distance. This is a position fac- 
tor of —1. We are left with two parameters that can be 
used to reduce aberrations, n and q. We will set the value 
of q to eliminate one of the aberrations and compare the 
remaining aberrations as a function of index. We will 
restrict our attention to moderate values of index. At 
large absolute values of index, the aberrations approach 
the same value independent of sign, but dielectric lenses 
with high index have significant reflection coefficients due 
to the impedance mismatch to free space. 

The usual ordering of the aberrations is from high- 
est to lowest in the order of r, the aperture coordinate. 
This is the ordering of most image degradation to least 
if one is forming images with significant lens aperture, 
but small to moderate image size, which is a common 
occurrence in applications. Thus, spherical aberration is 
an obvious target for elimination. However, there are no 
roots of C200 for values of index greater than one, which 
is why this aberration is referred to as spherical aberra- 
tion, since it appears to be inherent to spherical lenses. 
The usual practice is to eliminate coma (the next in line), 
and it so happens that the resulting lens has a value for 
the spherical aberration that is very near the minimum 
obtainable. Adjusting the shape factor, q, is often called 
lens bending. If we bend the lens for zero coma, that is 
find the roots of Cuo with respect to q we obtain 



9c 



(2n+l)(n-l) 
n+1 



(10) 



We plug this value for q and p = —1 into and plot the 
remaining three non-zero aberration coefficients as well as 
qc in Fig. |21 We note that there are two values of index 
where q — 1, which represent a plano-concave/convex 
lens. Setting H10|l equal to one we obtain. 



- n - 1 = 0. 



(11) 



the roots of which are the ubiquitous golden ratios, 
n = (f) ~ 1.62 and n = 1 — (/) ~ — 0.62,7J. We also note 
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FIG. 2: Top plot shows spherical aberration (black), astig- 
matism (green), field curvature (blue), and shape factor (light 
gray) as a function of index for a lens focusing an object at in- 
finity and bent for zero coma. Thin gray vertical lines indicate 
properties for lenses shown in ray tracing diagrams (bottom), 
meridional profile (left) and image spot (right). Incident angle 
is 0.2 radians and lenses are f/2. Index, shape factor, relative 
rms spot size, and spot diagram zoom are shown tabularly. 
In meridional profile, lens principle planes are shown as thin 
black vertical lines, and optic axis and Gaussian image plane 
are shown as blue lines. In spot diagram, Gaussian focus is 
at the center of blue cross hairs. 



that there is a window of index values near n = —0.7 
where both the spherical aberration and field curvature 
are small. There is no equivalent window in positive in- 
dex. 

Several ray tracing diagrams with both meridional rays 
and ray spot diagrams are shown for specific values of 
index in Fig. [5] The reference lens has index (f>, which 
is close to typical values used in visible optical lenses 
and near enough to n = 1 for reasonably low reflection. 
The lenses of negative index shown are in fact closer to 
n = — 1, which is the other index which permits perfect 
transmission, so this is a fair comparison. The negative 
index lenses all show significantly tighter foci than the 
positive index lens. 




FIG. 3: All as in Fig. |5| except the following. Lens is bent 
for zero spherial aberration. Coma is shown red. Solid and 
dashed lines indicate different solutions. Spot size, r^ms, is 
relative to bottom lens spot in Fig. |5| All spot diagrams are 
at the same scale. 



If we attempt to bend a lens with p — —1 to obtain 
zero spherical aberration we obtain the two solutions 

Is = ^ '—^ . 12 

n + 2 

These expressions have real values only for n < 1 /4, so an 
implementation of such a lens (embedded in free space) 
is not possible with normal materials. It is a surpris- 
ing and significant result that negative index permits an 
entire family of spherical aberration free spherical lenses 
that can focus a distant object to a real focus, Fig. O 
The solution with the negative sign in the expression for 
Qs (solid curves) has less coma for moderate negative val- 
ues of index, so ray tracing diagrams are shown for that 
solution. We note that at n = —1, the field curvature is 
also zero, thus this lens has only two of the five Seidel 
aberrations, coma and astigmatism. For a positive index 
reference we use the zero coma, n = cj) lens from above. 
Here again, negative index lenses achieve a tighter focus 
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FIG. 4: All as in Fig. |21 except the following. Lens configu- 
ration with object and image at finite positions and bent for 
zero spherial aberration and coma. Position factor is shown 
dark gray. Real image object pairs only occur when position 
factor is in shaded region, |p| < 1. Lens pairs are f/1.23, 
f/1.08, f/0.90 and have magnifications -1,-2,-3. In second to 
last spot diagram, horizontal (lOx) and vertical (lOOx) zoom 
are not equal. 



than a comparable positive index lens. 

Now we examine the case of |p| < 1, that is a real 
object and real image both at finite position. Since p and 
q are both free parameters, we can conceivably eliminate 
two aberrations. If we eliminate spherical aberration and 
coma the resulting lens is called aplanatic. It is a well 
known, though incorrect, result that a spherical lens can 
can only have virtual aplanatic focal pairs. The correct 
statement is that only negative index spherical lenses can 
have real aplanatic focal pairs. 

If we set C200 and Cuq to zero and solve for p and 
q, we obtain four solutions, the two non-trivial ones are 



given by 

n + 1 , , 

Psc = T 7, (13a) 

n — 1 

(7,c = ±(2n+l). (13b) 

We will focus on the solution with a minus sign for p and 
the plus sign for q. This solution has smaller aberrations 
for lens configurations that magnify an image. The other 
solution is better for image reduction. Inserting the ex- 
pressions (fT3^ into © we have plotted the two remaining 
non-zero coefficient as well as the values of Psc and qsc 
(Fig. 2J|. Ray diagrams are shown for lenses with mag- 
nifications of -1, -2 and -3. Also shown is a reference 
positive index lens for each. The reference lenses (which 
cannot be aplanatic) are of moderate index, 0, with the 
same magnification and f/# as the lenses they are com- 
pared to. They are bent for zero coma but also have 
spherical aberration near the minimum possible for the 
configuration. Again the negative index lenses produce 
superior foci. 

The lens of index —1 and magnification —1 is particu- 
larly interesting. At this index value the field curvature 
is also zero. This remarkable lens configuration has only 
one of the five Seidel aberrations, astigmatism. This is 
confirmed by ray tracing which shows a one dimensional 
"spot" at the image plane. This is perfect focusing in 
the sagittal plane. Perfect focusing also occurs in the 
meridional plane, in front of sagittal focus. 

One may ask why this asymmetric lens, g = — 1, per- 
forms so well in a symmetric configuration, p = 0. This 
lens can be equivalently viewed as a biconcave doublet 
with one component having index — 1 and the other hav- 
ing index 1, i.e. free space. Driven by this observation, 
we found that all biconcave doublets with arbitrary in- 
dices of ±n have identical focusing properties. The only 
observable difference is in the internal rays, which are 
always symmetric about the planer interface, but make 
more extreme angles at higher index magnitude. 

Fabrication of any of these negative index lenses is 
quite feasible using periodically structured artificial ma- 
terials. Current artificial material designs can operate at 
frequencies from megahertz through terahertzes], where 
there are numerous communication and imaging appli- 
cations. For example, lens antennas could benefit both 
by a reduction in aberrations, which translates directly 
into increased gain, and by a reduction of mass, afforded 
by low density artificial materials. Furthermore, these 
lenses are even easier to implement than a perfect lens, 
since they lack its severe structure period per wavelength 
requirements and are more tolerant to losses pj. Nega- 
tive index lenses at visible light frequencies may also be 
possible, by using photonic crystals, which have shown 
potential for negative refraction 10, ITj. 

Using the current optical system design paradigm, 
aberrations are minimized by combining elements with 
coefficients of opposite signjljl. However, more elements 
mean greater complexity and cost. Taking advantage 
of an expanded parameter space that includes negative 
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index can reduce the number of required elements — possibly even to one. 
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